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ABSTRACT 

We have analysed deep i?-band images, down to a limiting surface brightness of 
26.5 R-mag arcsec~ 2 (equivalent to ^28 B-mag arcsec~ 2 ), of 5 cD galaxies to determine 
the shape of the surface brightness profiles of their extended stellar envelopes. Both 
de Vaucouleurs i? 1 / 4 model and Sersic's i? 1 /™ model, on their own, provide a poor 
description of the surface brightness profiles of cD galaxies. This is due to the presence 
of outer stellar envelopes, thought to have accumulated over the merger history of the 
central cluster galaxy and also from the tidal stripping of galaxies at larger cluster 
radii. We therefore simultaneously fit two Sersic functions to measure the shape of 
the inner and outer components of the cD galaxies. We show that, for 3 out of our 
5 galaxies, the surface brightness profiles are best fit by an inner Sersic model, with 
indices n ~ 1 — 6, and an outer exponential component. For these systems, the galaxy- 
to-envelope size ratio is 0.1 - 0.4 and the contribution of the stellar envelope to the 
total i?-band light (i.e. galaxy + envelope) is around 60 to 80 per cent (based on 
extrapolation to a 300 kpc radius). The exceptions are NGC 6173, for which our 
surface brightness profile modelling is consistent with just a single component (i.e. no 
envelope) and NGC 4874 which appears to have an envelope with a de Vaucouleurs, 
rather than exponential, profile. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: formation - galaxies: 
fundamental parameters - galaxies: haloes - galaxies: structure 



1 INTRODUCTION 

First-ranked galaxies in clusters, also referred to as brightest 
cluster galaxies (BCGs) , are the brightest and most massive 
galaxies in the Universe. They are typically elliptical galax- 
ies (Lauer & Postman 1992). About 20% of BCGs appear 
to be surrounded by a large, low surface brightness envelope 
and are additionally referred to as cD galaxies (e.g. Dressier 
1984; Oegerle & Hill 2001). Such cDs reside only in clusters 
and groups, never in the field. Their existence and evolu- 
tion are intimately tied to the formation and evolution of 
the clusters themselves. The detection of these envelopes, 
however, is somewhat problematic. 

For a time, every elliptical galaxy was thought to have 
a stellar distribution whose projection on the plane of the 
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sky was described by de Vaucoulers (1948) R 1 ^ 4 law. This 
is reflected by the status of "law" that is ascribed to what 
is a highly useful, but nonetheless empirical " model" . How- 
ever, Lugger (1984) and Schombert (1986) have shown that 
all luminous, elliptical galaxies, including brightest cluster 
galaxies, have excess flux at large radii relative to their best- 
fitting i? 1 / 4 models. Moreover, today, it is known that only 
elliptical galaxies with Mb ~ —20.5 mag have E}/ A profiles 
(e.g. Kormendy & Djorgovski 1989). Brighter and fainter el- 
liptical galaxies are better described by Sersic's (1963) 7? 1//n 
model (see Graham & Driver 2005 for a review) with the in- 
dex n taking on values that are greater and smaller than 4, 
respectively (e.g. Caon, Capaccioli & D'Onofrio 1993; Young 
& Currie 1994; Graham et al. 1996; Graham & Guzman 
2003, and references therein). This then leads to the ques- 
tion as to whether the excess flux observed in cD galaxies 
(e.g. Feldmeier et al. 2002; Liu et al. 2005) is due to a distinct 
and separate halo of material, or is instead a manifestation 
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Table 1. Column 1: galaxy name. Column 2: host cluster name. Column 3: richness class from the catalogue of Abell (1958). Column 4: 
Type from Bautz & Morgan (1970). Column 5: galaxy redshift taken from the NASA/IPAC Extragalactic Database (NED). Column 6: 
physical scale, i.e. the distance in kpc that is equivalent to an angular distance of 1", calculated using a Hubble constant, Ho = 70 km 
s -1 Mpc - 1 . Columns 7 and 8: right ascension and declination of the galaxy. Columns 9 and 10: right ascension and declination of the 
nearby blank field. Column 11: worst seeing of each galaxy image. Column 12: radial extent of the data. 
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from the application of, or at least comparison with, an in- 
appropriate fitting function, namely the R 1 ^ 4 model. 

To address this question, and to derive both the size 
and flux ratio of any possible outer envelope relative to the 
inner galaxy component, obviously one should not apply the 
R 1 / 4 model to the inner light-profile. It would similarly be a 
mistake to simply fit an R}/ A model to any suspected outer 
halo. We have therefore set out to measure the shape of 
the projected stellar distribution through the simultaneous 
application of two R 1 ^" models to the light-profiles from 
deep exposures of galaxies reported to be cD galaxies. 

In this paper we analyse _R-band images for 5 cD galax- 
ies observed to a depth of ^jlr = 26.5 mag arcsec -2 , at the 3 
a level, with the main purpose of determining the shapes of 
the surface brightness profiles of their low surface brightness 
stellar envelopes. From this we also determine the galaxy- 
to-envelope size ratios and the envelope-to-total flux ratios 
within 300 kpc, and also when applying no outer trunca- 
tion to the best fitting modelfl These parameters can then 
be used to constrain models of cD galaxy, and host cluster, 
formation. This paper is arranged as follows. Section 2 de- 
scribes our observations and data reduction, including the 
ellipse fitting and surface brightness profile fitting method; 
Section 3 presents the results of the surface brightness pro- 
file fitting and discusses the best fit in each case. In Section 4 
we provide a brief review on the formation process of BCGs 
and then go on to discuss our findings. Finally, in Section 5 
we summarise our main results. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Data 

The light profiles presented in this paper are derived from 
broad i?-band images of 5 cD galaxies observed to a depth 
of hr = 26.5 mag arcsec -2 at the 3 a level. For elliptical 
galaxies, Worthey (1994) reports that B — R = 1.65 and 
Lauer & Postman (1994) find B — R — 1.51, and so our 
surface brightness limit is equivalent to a B-band depth of 
~ 28 mag arcsec -2 . This is ~3 mag arcsec -2 deeper than 
the study of BCGs presented by Graham et al. (1996). Our 



The radial extent of our data is typically ~ 100 kpc, and so we 
have to extrapolate to the truncation radius. 



galaxy sample was selected from two samples of BCGs: one 
observed by Lauer & Postman (1994), the other by Hill & 
Oegerle (1993). The selection criteria for these galaxies was 
that they (i) were classified as cD galaxies, (ii) had a limiting 
redshift z < 0.1, and (iii) were visible from La Palma in late 
April. Of the galaxies that met the selection criteria, we 
observed five; they are listed in Table 1. 

On each night at least four standard stars from the Lan- 
dolt (1992) list were observed, at varying airmass, in order 
to determine the photometric zero-point and the airmass 
extinction correction for the data. These have been applied 
to each image. Corrections for surface brightness dimming, 
galactic extinction and K-correction have also been applied. 
Throughout this paper, magnitudes are quoted using the 
Vega system. 

For each galaxy, one nearby blank field (see Table 1) 
was also observed for the purpose of estimating the sky- 
background free from any intracluster light. The observa- 
tions were taken in such a way that the object frames and sky 
frames were interleaved. The observations were made during 
dark time with the 1.0-m Jacobus Kapteyn Telescope (JKT) 
on the island of La Palma. These observations were taken on 
the nights of 2003 April 24-30 using the 2048x2048 SITe2 
CCD camera, which has 0.331 arcsec pixels and a field-of- 
view 11. 3x 11.3 arcmin 2 . The filter used was a standard Har- 
ris R band filter. The observations of both the galaxy fields 
and the blank fields consisted of 24 co-added exposures of 
900 seconds. The exceptions to this are NGC 3551 and NGC 
4874, where only 12 co-added exposures were made due to 
bad weather. 

Data reduction was performed within IRAF. All images 
were bias-subtracted, and then flat fielded using twilight 
flats. Images were then combined by degrading the best im- 
ages, by convolution with a Gaussian, to match the seeing 
conditions in the worst image, which was typically ~ 1" — 2" 
(see Table 1). There was no evidence for any fringing ef- 
fects in any of the images. The blank sky images were then 
used to estimate the sky-background and these values were 
subtracted from the galaxy images. In performing the sky- 
background estimation, ten areas of the blank sky image 
were used to calculate 10 medians. The sky value was then 
taken as the mean of these 10 medians, and the average devi- 
ation was adopted as the uncertainty in the sky background. 

Uncertainties in the estimation of the sky-background 
are the dominant source of error in determining the shape of 
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Figure 1. Differences in the sky background, Si, obtained from consecutive sky frames near each galaxy, and for the total dataset (bottom 
right panel). Each panel shows the associated sample standard deviation a = -\/) ]" ^/(n — 1), where n is 11, 22 or 88 depending on 
the panel. 



4 M. S. Seigar et al. 



Figure 2. Left: Images of the five galaxies in our sample. Right: Residual images derived by subtracting the azimuthally-symmetric 
(galaxy + envelope) model from the original images. 
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Figure 3. Surface brightness profile for GIN 478 with different analytic fits: (a) Sersic model; (b) double Sersic model; (c) Sersic + 
exponential model; (d) i? 1 / 4 + exponential model. Open circles at small radii are excluded from the fits. The root mean square (rms) 
scatter, A, is shown in the lower portion of each figure. The dashed lines indicate the extracted surface brightness profile with the sky 
uncertainty added and subtracted. We note that an arbitrary upper limit of fi c = 30 mag arcsec" 2 in our software is reached in panel 
(a), reflecting the inadequacy of a single i? 1//n function to describe the observed stellar distribution. 



the outer part of the galaxy surface brightness profiles. The 
above estimates of this error have been used in Section 2.3 to 
quantify the variation in the optimal parameters of the fit- 
ted models. Here we adopt an alternative approach to gauge 
the size of this error. Because of the way we interlaced our 
target and sky frames, the sky-background in each such pair 
of observations may differ. In Figure [1] we show histograms 
of the difference in the sky background between consecu- 
tive sky frames for each galaxy. Each offset in these his- 
tograms roughly represents twice the expected background 
offset between an individual target and sky frame, for which 
only half as much time has elapsed between their acquisi- 
tion. As can be seen, the median value for this difference 
is consistent with the sky frames having similar values for 
the background. The sample standard deviation, a, is only 
~0.012 mag arcsec -2 . A dominant monotonic drift in the 
background as one progresses through the night would pro- 
duce a non-zero median value in the distributions to either 
positive or negative values. Given N (=11 or 22) pairs of sky 
frames, the expected difference in the final sky-background 
level of each is given by a/yN. This amounts to 0.0026- 
0.0036 mag arcsec -2 uncertainty, consistent with the sky 
errors calculated using the method above. 



2.2 Extraction of surface brightness profiles 

We have employed and compared two methods for the ex- 
traction of the surface brightness profiles from the galaxy 
images. 

The first is isophotal ellipse fitting, performed using the 
Ellipse routine in IRAF which uses an iterative method de- 
scribed by Jedrzejewski (1987). Each isophote was fitted al- 
lowing for a variable position angle and ellipticity, but hold- 
ing the centre fixed. Foreground and background sources 
were masked out within the Ellipse routine. Given our in- 
terest in the outer light-profile, we chose to sample the sur- 
face brightness profile using a standard, albeit somewhat 
arbitrary, linear spacing. A logarithmic spacing would have 
generated a lot of data points, and hence more weight, at 
small radii. The Ellipse routine works in such a way, that 
ellipses are set up centered on the central cD galaxy. This de- 
fines several elliptical annuli. The surface brightness at each 
annulus is calculated as the average number of counts within 
that annulus. The uncertainty on the surface brightness is 
calculated as the rms between pixels divided by y/N — 1, 
where N is the number of pixels. The statistical error is 
therefore dependent upon the step size between annuli or 
ellipses. However, at large radii, the actual error on the com- 
puted surface brightness, is dominated by uncertainty in the 
sky background, and we discuss this later. 

The second method that we used determines the sur- 
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Figure 4. Same as Figure [3] but for NGC 3551. 
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Figure 6. Same as Figure [3] but for NGC 6173. The high value of H e ,2 in panel (b) is not a software limit. It is instead an indication 
that we may be merely fitting an inadequately subtracted sky-background with the outer R 1 '" function. 



face brightness profiles with the help of another isopho- 
tal fitting routine also written in IRAF (Jerjen, Kalnajs & 
Binggeli 2000; Barazza, Binggeli & Jerjen 2002). After fore- 
ground stars and neighbouring galaxies were removed from 
the image, a symmetrical 2-D model was reconstructed from 
the observed light distribution allowing isophotal elliptic- 
ity and position angle to vary with radius, but keeping the 
luminosity-weighted centre fixed. This fitting process was 
repeated iteratively until the residuals were minimised. The 
1-D surface brightness profiles were then calculated from 
the 2-D model by adopting mean values for ellipticity and 
position angles. 

Figure [2] shows the reduced images (left panel) and the 
residual images (right panel) of the 5 cD galaxies in our 
sample. The residual images were produced by subtracting 
an azimuthally symmetrical model from the original data. 

Both sets of surface brightness profiles, and the sub- 
sequent modeling, generated consistent results. From here 
on we refer to only the surface brightness profiles extracted 
using the IRAF task, Ellipse, presented in Figures [3] to [7] 
The dashed lines in Figures [3] to [7] represent the extracted 
surface brightness profile after adding (upper dashed line) 
and subtracting (lower dashed line) the uncertainty in the 
sky background. Note that this uncertainty was added to or 
subtracted from the image, and then the new surface bright- 
ness profile was derived. Because this uncertainty is applied 
to the image, the dashed lines sometimes do not necessarily 
fall either side of the data points. This is because it is pos- 
sible for the Ellipse routine to fit the image with slightly 



different ellipticities and position angles, and so a small ad- 
justment in the surface brightness results. 



2.3 Modelling the surface brightness profiles 

Our modelling of the surface brightness profiles employs the 
Sersic (1963, 1968) model for both the inner part of the 
galaxy and the outer stellar envelope. 

The Sersic R}/ n radial intensity profile can be written 

as 



I(R) = I e exp 



] -1 



f- 



(1) 



where I e is the intensity at the effective radius, R c , which 
encloses 50 per cent of the light. The factor b n is a function of 
the shape parameter, n, such that T(2n) = 2^(21%, b n ), where 
F is the gamma function and 7 is the incomplete gamma 
function (see Graham & Driver 2005). In the case where 
n — 1, the Sersic model is equivalent to an exponential, and 
when n = 4 it is equivalent to the J? 1//4 model. 

Initially, for all the galaxies, we attempt to fit the entire 
surface brightness profile with a single Sersic component. We 
also test the applicability of fitting an inner R 1/n model and 
an exponential model to the outer envelope (i.e., in a similar 
way to the case for disk galaxies, e.g., Andredakis et al. 1995; 
Seigar & James 1998; Graham 2001). We then go on to use 
a multitude of fits (not all shown) , which keep the inner and 
outer Sersic indices fixed at integer values between 1 and 4. 
Finally, we allow both of the Sersic indices to vary. 
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Figure 7. Same as Figure [3] but for UGC 9799. In panel (d), the outer exponential component does not contribute to the fit, and is 
thus not seen. The excess central flux seen in panels (b) and (c) are from the AGN. 



Corrections for the effects of seeing have been made us- 
ing the prescription given in Pritchet & Kline (1981). Due 
to the Gaussian nature of the JKT point spread function 
(psf) it is not necessary to consider more complicated seeing 
corrections. For any intrinsically radially symmetric inten- 
sity distribution, I(R), the observed seeing-convolved pro- 
file, I C (R), is 

poo 

I C (R) = cr-V* 2 / 2 " 2 / I{x)I {xR/a 2 )e~ x2 /2 " 2 xdx, (2) 
Jo 

where a is the dispersion of the Gaussian psf, which is 
equal to the full-width half maximum (FWHM) divided by 
a factor of 2.3548. The Iq term is the zeroth-order modified 
Bessel function of the first kind (e.g. Press et al. 1986). This 
approach to correcting light profile shapes for seeing was 
adopted by Andredakis, Peletier & Balcells (1995) and later 
by de Jong (1996). 

Due to the potential presence of partially depleted cores 
in luminous elliptical galaxies (e.g., Lauer et al. 1995; Gra- 
ham et al. 2003; Trujillo et al. 2004; Ferrarese et al. 2006 and 
references therein) , or instead the presence of multiple nuclei 
from semi-digested mergers, the innermost seeing-effected 
data points (~ 3", ~ 2 — 5 kpc) have been excluded from the 
fits. In comparison, the inner 10-20 kpc were excluded from 
the BCG analysis in Zibetti et al. (2005), where the FWHM 
was ~5 kpc. It is because of such features that one should 
not model integrated aperture magnitude profiles, in which 
every data point is effected/biased. Instead, one should fit 
the surface brightness profiles directly. Because of the AGN 



in UGC 9799, the inner 4 data points were excluded. Obvi- 
ously, given that we have excluded the most seeing affected 
data, the use of equation 2 to convolve our R 1/n models 
before fitting them to the observed light profiles is not so 
crucial. Deactivating the seeing correction has no significant 
affect. 

The best-fitting models were acquired using the sub- 
routine UNCMND from Kahaner, Moler & Nash (1989). At 
each iteration, the nonlinear Sersic functions are approxi- 
mated by a quadratic function derived from a Taylor series. 
The quadratic function is minimised to obtain a search direc- 
tion, and an approximate minimum of the nonlinear function 
along the search direction is found using a line search. The 
algorithm computes an approximation to the second deriva- 
tive matrix of the nonlinear function using quasi-Newton 
techniques. 

Common practice when fitting a model to some data 
set is to employ the use of the \ 2 statistic, such that the 
reduced-x 2 is given by, 



where m is the number of data points, Si is the ith residual 
(about the best-fitting model), Oi is the uncertainty on the 
ith data point and k is the number of parameters in the 
fitted model. 

Such an approach is highly desirable if one knows what 
the functional form of the underlying model is. However, 
when one does not know the form of the underlying model, 
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but instead has to assume some empirical function such 
as de Vaucouleurs model or Sersic's model, the use of the 
reduced-x 2 statistic can produce rather biased results. For 
example, at the centre of a galaxy, the signal-to-noise ratio 
is high, and thus the uncertainties (<Xj) that are assigned to 
the central data are small. These points thus have consid- 
erable weight in determining the best fit, obtained by min- 
imising the reduced-x 2 value. In the past, the reduced-x 2 
statistic has been used to fit Sersic bulges plus exponen- 
tial disks to spiral galaxy light profiles. However, due to the 
presence of additional (un-modelled) nuclear components, 
the bulge model and the simultaneously-fit disk model have 
been heavily biased (see Balcells et al. 2003). That is, be- 
cause the assumed model (bulge + disk) did not match the 
true underlying distribution (bulge + a disk + an additional 
nuclear component), the small uncertainties on the data at 
small radii heavily biased the fits to produce erroneous re- 
sults (see e.g., Schombert & Bothun 1987). 

We do not know the structural make-up of our 5 cD 
galaxies, they may contain a third component, such as a bar 
or a lens or indeed multiple nuclear components that we do 
not model. We therefore wish to avoid use of the x 2 statistic. 

Furthermore, another main concern is that we wish to 
quantify the stellar distribution of the suspected outer en- 
velope. In using equation [3j the (signal-to-noise)-weighted 
values of <Ji will act to erase the value or worth of the data 
at large radii. Moreover, in trying to gauge the influence of 
sky-background errors, likely to be a major source of uncer- 
tainty on the shape of the outer stellar distribution, the use 
of equation [3] would dilute the effect of adding and subtract- 
ing the uncertainty in the sky-background, and give one the 
false belief that their fitted models have less variance than 
they really ought. This is because correlated errors are not 
taken into account when computing the x 2 value. 

A common approach, which circumvents the above two 
problems, and which we have adopted, is to use the root 
mean square (rms) scatter 

A-^/SZ. (4, 

The results of our profile fitting are shown in Figures [3]- 
[3 The Sersic indices, effective radii and effective surface 
brightnesses are listed in Table 2. Uncertainties on the best- 
fitting parameters (Table 2) are obtained by repeating the 
fit to the surface brightness profiles after adding and sub- 
tracting the I a uncertainty in the sky-background level. It 
should be noted that the effective radii are model parameters 
that provide the optimal fit to the data, over the observed 
data range. They additionally reflect real, physical half-light 
radii only if the models can be extrapolated to infinity. 

While the introduction of additional free parameters in 
a fitted model will reduce the value of A, we show here that 
this is not the explanation for the improvement in the fit 
we obtain when changing from an i? 1//4 model to an R}^ n 
model. 

Our profiles have at least m — 125 measured points. 
R 1 / 4 + exponential fits have four free parameters, i.e. k = 4, 
and Sersic + Sersic fits have k — 6. For rando m residuals, 
the expected value of A scales as */ (m — k)/m. As a result, 
an improvement of only 0.8% would be expected for increas- 
ing the number of free parameters from 4 to 6. To have an 
improvement of 5% in A, a total of 16 free parameters are 



needed, and for a 10% improvement, 27 free parameters are 
needed. For three of our light profiles an improvement of 
more than 10% is shown when adopting a Sersic + Sersic fit 
over an R 1 ^ 4 + exponential fit. That is, the residuals about 
the i? 1 / 4 + exponential fit are not random, instead, there is 
structure indicating the inadequacy of this model and justi- 
fying the double Sersic model. NGC 3551 shows little, if any, 
improvement in its A, revealing that an R 1 ^" 1 + R 1 * 1 " 2 with 
ni = 4 and n% = 1 is appropriate for this galaxy. A sum- 
mary of the values of A found for each type of fit applied to 
each of our galaxies is shown in Table [3] We do not consider 
NGC 6173 here since only one component is necessary to 
model this galaxy (see Section 3). 

Unfortunately neither the ellipticity profile nor the po- 
sition angle profile yielded any clues to the transition from 
inner to outer component. Similarly, Zibetti et al. (2005, 
their Fig. 6) show that no change in the ellipticty profile is 
observed at the inflection of their surface brightness profile. 
In fact, there is no change in the behaviour of their elliptic- 
ity profile until ~160 kpc — a radius 8 times greater than 
their inner component's effective radius. 



3 RESULTS: GALAXY SURFACE 
BRIGHTNESS PROFILES 

GIN 478: 

The obvious structure in the residual profile of Figure 
[5^ reveals that this galaxy is not well described with a 
single-component R 1 ' n model. Figure [3p, in which two 
7? 1//n models have been fitted, shows that the outer part of 
this galaxy's surface brightness profile is almost exponential 
in nature. Fitting an outer exponential (Figure [3}j) shows 
that the profile shape of the best-fitting inner model does 
not change significantly, with a Sersic index n ~ 3. Both 
of these fits have residuals of A ~ 0.07 mag arcsec -2 . 
Fitting the inner part of the profile with an R 1 ^ 4 model 
and the outer part of profile with an exponential (Figure 
[3}1) increases the residuals by more than 10 per cent with 
A = 0.077 mag arcsec -2 . An R 1 / 3 + exponential model 
is the optimal fit when using integer values for the Sersic 
indices. 



NGC 3551: 

This galaxy is also not well modelled with a single R 1 ^" 
function (Figure [4^,), evidenced by the mismatch at small 
radii and the relatively large rms scatter A (c.f. Figure|4p, c 
and d). Figure 2p shows that the outer part of this galaxy's 
surface brightness profile is almost exponential in nature. 
Modelling the outer profile with an exponential shows 
that the inner R 1 ^" model's profile shape does not change 
significantly, with a Sersic parameter n ~ 6 (Figure |4j;). 
Both of these fits have residuals of A = 0.049 mag arcsec -2 . 
Figure 2Ji shows an R 1 ' 4 fit to the inner component and 
an exponential fit to the outer component, which also has 
a residual of A = 0.049 mag arcsec -2 A Sersic model with 
index n ~ 4 — 6 therefore appears to provide a good fit to 
the central galaxy in this case, with an outer exponential 
law again describing the envelope. The change from an 
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Table 2. Summary of the results of the double Sersic model fitting, except NGC 6173 for which we present the results from the single 
Sersic component fit. Column 1: galaxy name. Column 2: total absolute magnitude within a 300 kpc radius. Column 3 and 4: surface 
brightness )i e at the effective radius R c , respectively, from the inner Sersic fit. Column 5: inner Sersic index n\. Column 6: absolute 
magnitude within a 300 kpc radius for the inner component. Column 7 and 8: surface brightness at the effective radius, and this radius, 
for the outer Sersic component. Column 9: outer Sersic index n^. Column 10: absolute magnitude within a 300 kpc radius for the 
outer component. Given that the dominant source of error in the observed surface brightness profiles arises from uncertainties in the 
sky-background, all parameter errors are obtained by re-fitting the models to the images with the sky uncertainty added or subtracted. 



Galaxy 


Mtot 




Inner 








Outer 








< 300 kpc 




He,l 


m 


dinner 


Me,2 




112 


M outer 






(mag arcsec -2 ) 


(arcsec) 




< 300 kpc 


(mag arcsec -2 ) 


(arcsec) 




< 300 kpc 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


GIN 478 


-23.03±0.45 


24.08±0.13 


11.5±0.9 


3.0±0.1 


-21.63±0.03 


27.38±0.05 


151.0±49.1 


0.8±0.1 


-22.69±0.43 


NGC 3551 


-22.21±0.24 


26.27±0.49 


86.0±18.7 


6.4±0.7 


-21.50±0.03 


26.95±0.23 


209.2±49.8 


1.2±0.2 


-21.40±0.22 


NGC 4874 


-21.78±0.09 


21.44±0.23 


3.8±0.3 


1.1±0.2 


-18.50±0.02 


26.84±0.13 


230.0±26.1 


4.0±0.1 


-21.73±0.08 


NGC 6173 


-22.95±3.57 


30.00±5.47 


788±218 


14.0±7.5 












UGC 9799 


-22.26±0.45 


23.17±0.17 


15.7±3.1 


1.4±0.2 


-20.73±0.19 


25.00±0.19 


72.3±13.4 


1.0±0.3 


-21.95±0.29 



Table 3. Summary of the values of A found for each type of fit applied to 4 of our 5 galaxies. See Scction l2.2l 



Galaxy 


R l/n 


fll/4 + R l/4 


A (mag 
H l/3 + R l/3 


arcsec 2 ) 
i? 1 / 4 +Exp 


ftl/n_|_Exp 


R l/n + R l/n 


GIN 478 


0.277 


0.118 


0.076 


0.077 


0.069 


0.068 


NGC 3551 


0.080 


0.050 


0.050 


0.049 


0.049 


0.049 


NGC 4874 


0.053 


0.042 


0.046 


0.057 


0.054 


0.040 


UGC 9799 


0.097 


0.130 


0.105 


0.130 


0.080 


0.080 



inner Sersic parameter of n = 4 to n = 6.4 results in an 
increase in R e from 37 to 86 arcsec, i.e. more than a factor 
of two. This is the only galaxy for which this degeneracy is 
seen, and this is not typical of our galaxies in general. 



NGC 4874: 

The residual profile in Figure [5^, resembles that seen in 
Figures [5^ and [3K revealing additional structure that a 
single-component Sersic model cannot describe. Curiously, 
Figure [5p shows that fitting the system with a double 
Sersic model provides an outer Sersic index of n ~ 4, which 
is not seen in any of the other galaxies. Moreover, the 
inner galaxy seems well characterized by a profile that is 
close to exponential. If we force the outer model to have 
an exponential profile (Figures [5J; and d), the residuals 
increase significantly from A = 0.040 mag arcsec -2 to 
A =0.054-0.057 mag arcsec -2 . The inner component of 
NGC 4874 has a notably small effective radius of ~4 kpc. 
This is perhaps not unusual as Gonzalez et al. (2005) also 
have a number of BCGs with -Runner < 5 kpc. 



NGC 6173: 

The residual profile in Figure [6^, demonstrates that a 
single-component fit may be adequate for this galaxy, 
although a large Sersic index (n ~ 14) is required, and the 
values of /i c and R c are unreasonably large. Indeed, the 
value of /i c = 30 mag arcsec -2 is an artificial upper limit 
used in the fitting code. Including an outer component 
decreases the residual from A = 0.077 to A = 0.063 (Figure 



[6p and [6fc) and requires an inner Sersic index of n ~ 7. 
However, given the near constant surface brightness of this 
outer component over the radial range for which we have 
data, at first glance it looks as if we failed to adequately 
subtract the sky background. The additional subtraction 
of the 1 a uncertainty in the sky background results in 
the disappearance of this near constant surface brightness 
component seen in Figures [6p, c and d, and the optimal 
fit is actually a 1-component fit, with n = 8.7, R c = 141.2 
arcsec and jj, c = 25.1 mag arcsec -2 . We thus do not claim to 
have detected a distinct, outer halo in this systerrQ, which 
may therefore not be a cD galaxy. That is, this system may 
simply have been mis-classified as a cD galaxy because of 
excess flux at large radii compared to an R 1 ^ 4 model. 

UGC 9799: 

Figure [7p shows the double Sersic fit to this galaxy. The 
inner and outer components are both modelled with n ~ 1, 
suggesting that this galaxy's surface brightness profile may 
be well characterized by a double exponential. Keeping the 
outer component modelled with n — 1 (i.e. exponential) and 
the inner component modelled with a free value of n does 
not significantly change the residuals (Figure Et) and shows 
that the inner component still has a Sersic index n ~ 1. 
Not surprisingly, the affect of fitting this galaxy with a dou- 

2 The lack of a need for an outer component is not unprecedented. 
Gonzalez et al. (2000) present the light-profile for the BCG+halo 
in Abell 1651, and although it extends to an impressive 670 ft. -1 
kpc, there is no sign of a transition radius and they show it is 
remarkably well fit by a sing le R 1 / 4 model. 
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Figure 8. HST/WFPC2 F814W images of the inner 5" of UGC 9799. In all images north is up and east is left. Top left: A high contrast 
image showing signs of fine structure to the north-east of the nucleus of UGC 9799. Top right: A low contrast image showing the AGN. 
Bottom left: An unsharp-masked image, highlighting the fine structure to the north-east of the nucleus. Bottom right: A structure map 
(see Pogge & Martini 2002) highlighting the AGN. 



ble exponential (not shown) does not change the residuals 
significantly. It is unusual though to see such a large ellipti- 
cal galaxy with a low value for n, since luminous early-type 
galaxies usually have high Sersic indices (Caon, Capaccioli 
& D'Onofrio 1993; Graham et al. 1996). 

FigureOshows a series of /-band (F814W) HST/WFPC2 im- 
ages of UGC 9799 (taken from the HST archive and observed 
as part of program number SNAP-8683 PI: van der Marel). 
Images are shown with different stretches, highlighting both 
low-surface brightness features (left), and the fact that this 
galaxy has an active nucleus (right), which is classified as 
Seyfert 2 in the NASA Extragalactic Database (NED). The 
image in the bottom left of Figure [8] is a type of residual 
image, created by fitting the smooth galactic starlight with 
elliptical isophotes and subtracting the original image from 
this model. Fine structure emerges using this technique. In 



this case low surface brightness features to the north-east of 
the nucleus have been revealed. The image in the bottom 
right of Figure [8] is a structure map (of the kind presented 
by Pogge & Martini 2002) and again highlights the AGN. 

For our ground-based surface brightness profiles, the 
inner few points (determined from the size of the relevant 
seeing disk) are ignored when performing the analytical fits. 
As a result, the bulk of the AGN contribution will also be 
ignored, as this will be the same size as the seeing disk. The 
fine structure seen in Figure 7 is also common in central clus- 
ter galaxies and can be interpreted as features, which appear 
as a result of the merger processes involved in the formation 
of cD galaxies. From a morphological point of view, UGC 
9799 shows nothing that would not be expected for a cD 
galaxy, apart from an inner exponential fit. Although, this 
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was also observed in NGC 4874 and thus may be more com- 
mon in cD galaxies than previously thought. 



4 DISCUSSION 

Constraining the surface brightness profiles of the faint en- 
velopes in cD galaxies is important for constraining cur- 
rent models of cD galaxy formation and cluster dynam- 
ics. There is evidence from the globular cluster population, 
and the near-infrared galaxy luminosity function, to sug- 
gest that BCGs experienced their mergers long ago (Jordan 
et al. 2004; Ellis & Jones 2004), yet the presence of (un- 
erased) tidal streams (e.g. Gregg & West 1998; Trentham & 
Mobasher 1998) would appear to favour a more recent for- 
mation epoch. One of the earliest studies of the low surface 
brightness haloes of cD galaxies was performed by Carter 
(1977), who found that the total luminosity did not con- 
verge even at a radius of 300 kpc. This was later confirmed 
by Lugger (1984), who additionally found that a de Vau- 
couleurs R 1 ^ 4 model consistently underestimated the sur- 
face brightness profiles of cD galaxies at large radii (see also 
Schombert 1986). These observations have been confirmed in 
recent years by deep imaging of central cluster galaxies (e.g. 
Feldmeier et al. 2002, 2004a, b; Lin & Mohr 2004; Adami 
et al. 2005; Gonzalez, Zabludoff & Zaritsky 2005; Kemp et 
al. 2005; Krick, Bernstein & Pimbblet 2006; Liu et al. 2005; 
Mihos et al. 2005; Zibetti et al. 2005) which clearly indicate 
the presence of an extended stellar envelope, albeit relative 
to an inner R 1 ^ 4 component. The need for two Sersic com- 
ponents, rather than one, is clearly illustrated in Gonzalez, 
Zabludoff & Zaritsky (2003, their Figs. 1 & 2). Furthermore, 
it has also been shown that the Petrosian properties of cen- 
tral cluster galaxies display distinct properties, which can 
be interpreted as an indicator of cD galaxy haloes, indepen- 
dent of an assumed R 1 ^ 4 light-profile for the central galaxy 
(Brough et al. 2005; Patel et al. 2006). 

We note that our photometry alone does not indicate 
whether the outer 'component' is physically distinguished 
from the inner 'component'. In general, however, a single 
Sersic function provides a good fit to ordinary elliptical 
galaxies, with little or no structure in their resultant residual 
profiles. The structure seen in the residual profiles for 4 of 
our 5 cD galaxies is reminiscent of that seen when fitting a 
single Sersic function to a spiral galaxy, and is suggestive of 
two distinct components. The failure of a single Sersic func- 
tion to match the observed stellar distribution is physical 
evidence that four of our objects are different from ordinary 
elliptical galaxies, but we caution that they may still be sin- 
gle physical entities. 

A small fraction of the ICL, and thus envelopes around 
cD galaxies, probably originates from stars which have been 
gravitationally ejected by supermassive black hole binaries 
at the centres of elliptical galaxies within the cluster (e.g. 
Holly-Bockelmann et al. 2006) . Recently, Graham (2004) has 
shown that the central stellar mass deficit in "core" galaxies 
— thought to have formed in "dry" mergers — is roughly 
0.1 per cent of their total stellar mass (see also Ferrarese 
et al. 2006). This mass deficit is roughly equal to the (com- 
bined) mass of the central black hole in ellipticals, and is 
also consistent with theoretical predictions on the orbital 
decay of binary black holes (Ebisuzaki, Makino & Okumura 



1991; Milosavljevfc & Merritt 2001). Most recently, in high- 
precision, iV-body simulations, Merritt (2006) has shown 
that virtually all of the mass deficit is generated during the 
rapid, initial phases of binary formation, not after the binary 
becomes hard. He obtains mass deficits on the order of the 
mass of the binary's larger black hole, and so one can ex- 
pect the intracluster light (from this mechanism) to roughly 
equal ~ 0.1 per cent of the cluster light in spheroid^. 

Observations of multiple nuclei within brightest clus- 
ter galaxies (e.g. Hoessel & Schneider 1985; Postman & 
Lauer 1995; Seigar, Lynam & Chorney 2003) and low surface 
brightness tidal features (e.g. van Dokkum 2005) are consid- 
ered strong evidence that massive galaxies are growing at the 
centres of rich clusters by accreting their less massive neigh- 
bours (Hausman & Ostriker 1978). Such merger events are 
also thought to be partly responsible for the formation of 
extended (low surface brightness) envelopes and intraclus- 
ter light (Ostriker & Tremaine 1975; Ostriker & Hausman 
1977; Hill & Oegerle 1998; Moore et al. 1996; Muccione & 
Ciotti 2004; Willman et al. 2004). This is demonstrated in 
recent semi-analytical models (e.g., Purcell, Bullock & Zent- 
ner 2007) and N-body simulations (e.g., Conroy, Wechsler & 
Kravtsov 2007). Moreover, close galaxy-galaxy encounters 
can strip stars from deep within their potential well. These 
stars may then be liberated by the overall cluster tidal field 
(Merritt 1984; Moore et al. 1996), to become what is known 
as the intracluster light (ICL: Zwicky 1951; Welch & Sastry 
1971; Oemler 1973; Thuan & Kormendy 1977). Furthermore, 
the large numbers of ultra compact dwarf galaxies found in 
galaxy clusters suggests that this process may be rather ef- 
ficient (e.g., Bekki et al. 2003; Cote 2005; Drinkwater et al. 
2005; Gnedin 2003; Mieske, Hilker & Infante 2005). It is 
also likely responsible for the existence of intracluster: plan- 
etary nebulae (Arnaboldi et al. 2004; Aguerri et al. 2005a; 
Feldmeier et al. 2004b; Gerhard et al. 2005); red giant stars 
(Ferguson, Tanvir & von Hippel 1998; Durrell et al. 2002); 
novae (Neill, Shara & Oegerle 2005); and supernova (Gal- 
Yam et al. 2003). 

Our results reveal, given the validity of these processes, 
that they result (in three instances) in an exponential-like 
distribution of stars around the central dominant galaxy 
More precisely, for three of our five galaxies, an inner Sersic 
model plus an outer exponential model provides a good fit 
to the data. In one galaxy, NGC 6173, no outer exponen- 
tial model is required, and in NGC 4874, the outer light 
profile is best described with an R 1 ^ 4 model rather than an 
exponential model. 

If the envelopes associated with cD galaxies trace a sur- 
rounding dark matter halo, then one might expect them to 
be described by a Sersic function with n around 2.5 to 3. 
This is because hierarchical ACDM simulations produce a 
near universal profile shape for dark matter halos on all 
scales, the projection of which is well described by a Sersic 
R 1/n model with n ~ 3 (Lokas & Mamon 2001; Merritt et 
al. 2005, 2006). 

Curiously, iV-body simulations of cold collapses (and 
disk galaxy mergers) also result in haloes (and merger rem- 



3 This is an upper estimate because in gas rich mergers, gas fa- 
cilitates the decay of the binary black hole, and consequently less 
stars are ejected. 
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nants) having, in projection, i? 1//3 -like profiles (e.g., Will- 
man et al. 2004, their Fig. 7; Aceves, Velaquez & Cruz 2006; 
Merritt et al. 2006). Interestingly, a closer inspection of old 
data (e.g. Figures 4-6 in van Albada 1982) reveals obvious 
and systematic deviations from the R 1 ^ 4 model in the sense 
that an R 1 ^" model with n < 4 provides a better fit to the 
data presented there. If the halo or intracluster light around 
cD galaxies traces the dominant dark matter potential, one 
does not expect this envelope to be described with an R 1 ^ 4 
model. Of course, baryonic processes may well result in a 
different stellar distribution to the dark matter distribution, 
and studies suggest this can lead to a flattening of the in- 
ner density profile (Nipoti et al. 2004), which is effectively 
equivalent to a reduction in the Sersic index. 

Demarco et al. (2003) have analysed the distribution 
of X-ray gas in a sample of 24 real galaxy clusters. They 
found it was well described by a Sersic function having val- 
ues 0.8 < n < 2.3. With the exception of NGC 4874, this 
range is in good agreement with the values reported here for 
the outer component of our cD galaxies, and suggests that 
this envelope is indeed tracing the (azimuthally-averaged) 
intracluster light. Moreover, our exponential-like outer pro- 
files match the exponential ICL profile in Abell 3888 (Krick 
et al. 2006). When more observations become available, it 
will be interesting to see whether the distribution of the in- 
tracluster stellar probes, such as planetary nebula and glob- 
ular clusters, follow an exponential or an _R 1//4 radial distri- 
bution. It will also be interesting to know if the intragroup 
light (e.g., Da Rocha & de Oliveira 2005; Faltenbacher & 
Mathews 2005; White et al. 2003) behaves in a similar or 
different manner. 

The mean BCG+ICL light-profile obtained from the 
stacked cluster image reported in Zibetti et al. (2005) is plot- 
ted using a linear radial axis in Zibetti & White (2004, their 
Fig. lfl One can immediately see by eyeball examination 
that the outer light-profile is well approximated by an expo- 
nential (i.e., a straight line in that figure). This is in good 
agreement with our independent data and the (cluster halo) 
X-ray data from Demarco et al. (2003), but at odds with 
the J? 1 / 4 model used in Gonzalez et al. (2005), and at odds 
with a projected NFW model. We do however note that the 
shape of the ICL profile in Zibetti et al. (2005) is different, 
suggestive of a Sersic index greater than 1. This difference 
arose from their new corrections for mask incompleteness 
and their new method of determining the sky background. 
The latter involved simultaneously fitting an NFW model 
for the ICL and some constant value for the sky background 
level. The problem with such an approach is that the fitted 
constant effectively modifies the real ICL profile to produce 
(as best as it can) an NFW profile, even when the real ICL 
may not have such a form. Fitting an i? 1//4 model from 150 
to 500 kpc, Zibetti et al. (2005, their section 5.1) report an 
effective radius _R C of 250-300 kpc for the ICL, somewhat 
larger than the values we observe in our sample (Table 4). 



4 Zibetti & White (2004) use 654 clusters, while Zibetti et al. 
(2005) use 683. 



4.1 Relative contribution of the stellar envelope 
to the total luminosity of cD galaxies 

We use the analytical fits from Table [2] to determine the 
relative contribution of the stellar envelope to the total lu- 
minosity of the galaxy plus envelope, i.e. the envelope-to- 
total ratio. We compute two estimates of this flux ratio. 
One of these, ( -E/T) 300, comes from truncating the models 
at a radius of 300 kpc, as done in Gonzalez et al. (2005). 
Our second estimate assumes no truncation, and is denoted 
(E/T) to t- Due to the extended nature of the envelope, our 
(E/T) to t values are somewhat different when compared to 
our (E/T)soo values, although only by ~10 per cent at most. 
Both quantities are listed in Table 4. It should be noted that 
since the radial extent of our surface brightness profiles is 
~100 kpc, although 223 kpc for GIN 478, that we have to 
extrapolate to estimate both the (E/T)soo and (E/T) to t val- 
ues. 

In most cases the extended stellar envelope contributes 
around 60 to 80 per cent of the total i?-band luminosity 
when no truncation radius is applied and the models are 
extrapolated to infinity. When the profiles are truncated at 
300 kpc, the envelope contributes between ~45 to 80 per 
cent. The exception is NGC 4874 (the only galaxy with an 
R 1/4 envelope) which has E/T ~ 95-98 per cent. These 
flux ratios lie in the same range found by Gonzalez et al. 
(2005), who reported E/T ratios of around 0.9 but as low 
as 0.4 (their Figure 7). However, Zibetti et al. (2005) find 
a slightly lower value. They find that the ICL contributes 
10.9 percent to the total cluster light and the central galaxy 
contributes 21.9 per cent. This is equivalent to an E/T ratio 
of ~33 per cent. However, taking into account the errors in 
the measurements, E/T ratios in the range 25 to 45 per 
cent are allowed. The high end of their range is therefore 
consistent with the results found here. 

For the three cDs best described with an R 1 ^" galaxy 
plus exponential envelope, the galaxy-to-envelope size ratio 
(given by the ratio of the effective radii Re,i/Re,2) ranges 
from ~0.1 to ~0.4. In contrast, using double R 1 ' 4 models, 
Gonzalez et al. (2005) report ratios of ~0.1 down to ~0.025, 
i.e. envelopes 10 to 40 times larger in size than the inner 
component. One of our remaining two galaxies appears to 
have no distinct envelope, and the other has an R 1 ^ 4 enve- 
lope 60 times greater in size than the central galaxy. 



5 SUMMARY 

We have observed 5 cD galaxies to a depth of /j,r = 26.5 mag 
arcsec -2 , and we have determined the shapes of the surface 
brightness profile of their outer stellar envelopes. 

The results of previous attempts to model the intraclus- 
ter light or extended stellar envelope suggested that a uni- 
versal model applies. For example, Gonzalez et al. (2005) re- 
port that both the central part of the cD and the intracluster 
light are both well described by an R 1 ^ 4 surface brightness 
model. In general, previous studies of this kind have only 
tried fitting a surface brightness model of one kind. Our 
approach of fitting a Sersic model to the extended halo pro- 
vides a means to actually measure, rather than pre-ordain 
the actual stellar distribution, albeit within the confines of 
the Sersic model. 
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Table 4. Results using the double Sersic model to derive ratios of physical parameters. Column 1: galaxy name. Column 2: effective 
radius of the central galaxy. Column 3: effective radius of the 1CL or envelope. Column 4: galaxy-to-envelope size ratio, given by the ratio 
of the effective radii R e ,l/ Re,2- Column 5 and 6: envelope-to-total ratio calculated by extrapolating the profiles to infinity, (E/T)tot, 
and 300 kpc, (E/T) 300 . 



Galaxy 


Re,l 


Rc.2 


Re,l/Re,2 


(E/T) tot 


(E/T) 300 




(kpc) 


(kpc) 








1 


2 


3 


4 


5 


6 


GIN 478 


19.9±1.6 


261.2±84.9 


0.076±0.019 


0.82±0.13 


0.73±0.13 


NGC 355f 


53.3±11.6 


129.7±30.9 


0.4fl±0.072 


0.59±0.04 


0.47±0.05 


NGC 4874 


1.8±0.1 


108.1±12.3 


0.0f7±0.002 


0.98±0.01 


0.95±0.01 


NGC 6173 


1073.1±296.8 










UGC 9799 


10.5±2.0 


48.4±9.0 


0.2f7±0.042 


0.77±0.14 


0.76±0.14 



Our analysis suggests that the surface brightness pro- 
files of cD galaxies (including their envelopes) are best mod- 
elled by a double Sersic function. While an inner R 1/A model 
is sufficient for some cDs, we have found that the inner Sersic 
index can vary significantly from object to object (from 
n ~ 1 to n ~ 7). An outer exponential model seems ap- 
propriate for three of our four, 2-component systems. One 
galaxy (NGC 4874) appears to have an R 1 ^ envelope. Typ- 
ically, when present, the envelope contributes around 60 to 
90 per cent of the total (galaxy + ICL) light, and the galaxy- 
to-envelope size ratio is ~0.1 to ~0.4. 
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